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Abstract 

The multiblock three-dimensional Navier-Stokes method 
PAB3D was utilized by the Component Integration Branch 
(formerly Propulsion Aerodynamics Branch) at the NASA- 
Langley Research Center in an international study sponsored 
by AGARD Working Group #17 for the assessment of the 
state-of-the-art of propulsion-airframe integration testing 
techniques and CFD prediction technologies. Three test 
geometries from ONERA involving fundamental flow physics 
and four geometries from NASA-LaRC involving realistic 
flow interactions of wing, body, tail, and jet plumes were 
chosen by the Working Group. An overview of results on four 
(1 ONERA and 3 LaRC) of the seven test cases is presented. 
External static pressures, integrated pressure drag and total 
drag were calculated for the Langley test cases and jet 
plume velocity profiles and turbulent viscous stresses were 
calculated for the ONERA test case. Only selected data from 
these calculations are presented in this paper. The complete 
data sets calculated by the participants will be presented 
in an AGARD summary report. Predicted surface static 
pressures compared favorably with experimental data for 
the Langley geometries. Predicted afterbody drag compared 
well with experiment. Predicted nozzle drag was typically 
low due to over- compression of the flow near the trailing 
edge. Total drag was typically high. Predicted jet plume 
quantities on the ONERA case compared generally well with 
data. 

Introduction 

The AGARD Fluid Dynamics Panel sponsored Work- 
ing Group (WG) #17 for the investigation of the aero- 
dynamics of 3-D aircraft afterbodies. This group extended 
the work of Working Group #8 which reported on the aero- 
dynamics of 2-D afterbodies in 1986, ref. 1. The panel rec- 
ognized the extremely complex nature of the flow in the aft 
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region of fighter aircraft due to the various interactions 
of the flows about the fuselage, empennage, and nozzle 
boattails. These flows can involve such features as separated 
flow, shock/boundary layer interactions and multiple jet 
exhaust plumes. The panel, in addition to seeking to develop 
recommendations for improvements in propulsion testing 
techniques, sought to continue the evaluation of state-of- 
the-art computational methods for their ability in predicting 
these complex afterbody flows. 

The Working Group members selected seven test geome- 
tries for the evaluation program. A combination of single 
and twin-jet configurations were selected emphasizing dif- 
ferent aspects of afterbody flows. Three test configurations 
proposed by ONERA (designated A.x) were chosen for the 
evaluation of prediction techniques for basic flow physics. 
These configurations had extensive flow measurement data- 
bases such as off-body flow quantity surveys with both hot- 
and cold-jet simulations (A.2). Four Langley test cases 
(designated B.x) were chosen for being more realistic aero- 
dynamic configurations that would involve complex attached 
and separated external flows interacting with simulated jet 
plumes. The data acquired on these models included surface 
static pressure distributions and force balance data on the 
model afterbody. 

Several issues impair accurate code calibration/ validation 
efforts. Excepting the validity of the computational fluid dy- 
namics (CFD) itself, knowledge of the actual model shape, 
accuracy of the instrumentation, interference from facility 
related factors, and the accuracy of setting the facility flow 
parameters all contribute to uncertainties in the compar- 
isons. Most of the test geometries for the WG#17 study 
were well documented. Several of the models were measured 
after the tests to determine the manufactured shape. All 
the test cases had well defined free-stream and jet operating 
conditions. 

Highlights of results of the Component Integration Branch 
(CIB), NASA-LaRC (formerly Propulsion Aerodynamics 
Branch) involvement in the WG#17 program are presented. 
The multiblock three-dimensional Navier-Stokes method 
PAB3D was used for all the calculations. Jet plume velocity 
profiles and turbulent viscous stresses were calculated for the 
A.2 ONERA test case. External static pressures, integrated 
pressure drag and total drag were calculated for the Langley 
test cases B.l, B.2 and B.3. Calculations were made on B.4 
by CIB but are not presented here. 
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Nomenclature 

maximum body cross-sectional area, 

0.0273 m 2 

wing reference area, 0.4286 m 2 
drag coefficient, goo A ma x (BA,B.2); 

cfcj (B 3) 

pressure coefficient, ^~£°° 

body maximum diameter (ONERA A. 2), m 

axial force along body axis, N 

length of body (0.18212, B.1,B.2; 

0. 17474, B.3),m 

free-stream Mach number 
nozzle pressure ratio, ^ 
static pressure, Pa 
jet total pressure, Pa 
free-stream static pressure, Pa 
free-stream dynamic pressure, Pa 
Reynolds number 

jet total temperature, (typ. cold jet 
300 deg), K 

free-stream velocity, m/sec 
axial velocity, m/sec 

Reynolds stress velocity components, m/sec 
axial distance, m 

lateral distance from model centerline, m 

vertical or radial distance from model 
centerline, m 

flow angle- of- at tack, deg 

angular location of pressure orifices, deg 


Subscripts 

A afterbody component contribution 

CL centerline 

N nozzle component contribution 

P pressure force contribution 

T total forces from components 

t total conditions 

HT horizontal tail component contribution 

VT vertical tail component contribution 


Computational Procedure 
Governing Equations 

The code used by CIB at Langley was the general 3-D 
Navier-Stokes method PAB3D. This code has several com- 
putational schemes and different turbulence models that can 
be utilized, as described in more detail in references 2 and 3. 
The governing equations are the Reynolds- averaged simpli- 
fied Navier-Stokes equations (RANS) obtained by neglect- 


ing all streamwise derivatives of the viscous terms. The re- 
sulting equations are written in generalized coordinates and 
conservation form. The implementation of the full three- 
dimensional viscous stresses are reduced to thin-layer vis- 
cous assumptions, although full Navier-Stokes simulation is 
an option. The diffusion terms are centrally differenced and 
the inviscid flux terms are upwind differenced. Two finite 
volume flux-splitting schemes are used to construct the con- 
vective flux terms. The Roe upwind scheme with third order 
accuracy is used in evaluating the explicit part of the govern- 
ing equations and the van Leer scheme is used to construct 
the implicit operator. 

The user will typically utilize the Roe scheme procedure 
to sweep streamwise through the computational domain and 
the van Leer scheme for the solution of the cross-plane (i.e., 
f=constant) of a three dimensional problem. A single-cell 
wide, two-dimensional mesh defined with the i direction of 
the grid oriented in the conventional streamwise direction 
will converge more slowly using the Roe relaxation solution 
scheme compared to solving the equivalent problem with the 
van Leer scheme. Therefore the i and j directions of the 
2-D mesh are swapped allowing the entire flowfield to be 
solved implicitly with each iteration. The explicit sweep is 
not used since only one cell exists in the i direction. The 
implicit scheme has the potential of a much higher rate of 
convergence and typically provides a solution in less time. 

Several near-wall models and compressibility corrections 
are available to be used with several formulations of linear 
two-equation k-e turbulence model, (e.g. Standard, Jones & 
Launder, Yang & Shih). The k-e turbulence model equa- 
tions are uncoupled from the RANS equations and can be 
solved with a different time step than that of the princi- 
ple flow solution. The latest version of the code (V.13) has 
the capability to provide Algebraic Reynold’s Stress turbu- 
lence simulations, though this was not implemented in the 
code at the time of this study. All flow solutions were de- 
veloped with PAB3D-V10 using the Jones & Launder high 
Reynolds number formulation of k-e. The jet plume turbu- 
lence simulation also utilized the high-Reynolds number for- 
mulation with the compressible dissipation of Sarkar, ref. 4. 
The damping function of Launder and Sharma was used to 
control the near-wall behavior of k and e was set to zero at 
solid surfaces, ref. 5. 

Solution Process 

The calculation of laminar flow solutions require the 
memory usage of 15 words per grid point. Turbulent flow 
solutions using the linear two-equation k-e model requires 
20 words per grid point. The code speed was 43 //sec per grid 
point (Cray 2 time) for solving turbulent flow simulations. 

Several parameters were used to gauge solution conver- 
gence for the A. 2.2 case. Afterbody drag and nozzle mass 
flow, calculated by the performance package, ref. 6, and so- 
lution residual were used to determine the solution status at 
the coarse, medium, and fine grid levels. Typically mass flow 
variance of less than 0.1 percent and afterbody drag variance 
of less than 0.25 percent for several hundred iterations were 
achieved. Similar indicators for convergence were used for 
the B.l cases. External surface static pressure distributions 
were monitored for an indication of convergence for the B.2 
and B.3 cases. 
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Grid Topology and Block Interfaces 

All grids utilized H-H or H-0 type mesh topologies and 
had block dimensions that were multiples of 4. The ONERA 
A.2 case was gridded with a single cell wide 5 degree wedge 
grid with the streamwise flow direction oriented along the j 
index to utilize the implicit flow solver in the code for faster 
solution convergence. The B.l mesh had the same topology 
as the A.2 mesh. The conventional grid orientation was used 
for the Langley B.2 and B.3 cases, i.e. i index, along which 
the solution relaxation scheme occurs, is oriented in the 
streamwise direction. Block faces were matched one-to-one 
or integer-to-one, except for the A.2 case where some block 
face cell mis-matches occurred for gridding convenience. 

The conservative patch interface package of Pao & Abdol- 
Hamid, ref. 7, enables the code to properly transmit infor- 
mation between mis-matched block interfaces. Integer-to- 
one interfaces are considered a subset of the arbitrary block 
interface and do not need to be specified as such to the 
patching code. The patching program writes a connectiv- 
ity database as a preprocessor prior to execution of the flow 
solver. Each entry to the patch database contains cell face 
areas and indices relating that cell with all other cells that 
will share momentum flux information. The database in- 
formation is automatically re- allocated internal to the code 
during mesh sequencing. As a result, each block can be se- 
quenced at different levels and the correct interface informa- 
tion is maintained at the cell level. However, it is important 
to note that features in the flow developed on one side of an 
interface should not be obliterated on the other side due to 
too severe a grid density mis-match. 

Third-order continuity in transmitting the fluxes across 
block boundaries is maintained by the code; lower order 
continuity may be specified by the user if required. Equal 
cell size spacing on either side of an interface in directions 
normal to the interface must still be maintained regardless 
of the mesh sequencing level of the block. 

Boundary Conditions 

For this study, solid walls were treated as no-slip adiabatic 
surfaces. The solid wall boundary condition was satisfied by 
setting the momentum flux of the solid wall cell face to zero. 
The boundary conditions used for the internal nozzle flow 
path were fixed total pressure, total temperature and flow 
angle at the plenum block inflow face. A boundary con- 
dition for the Riemann invariants along the characteristics 
was specified for the external freestream inflow face and the 
lateral freestream outer boundary of the flow domain. An 
extrapolation boundary condition was applied on the down- 
stream outflow face where both the free stream and the noz- 
zle plume exit the computational domain. The axisymmetric 
flow assumption for the single-cell grids was implemented by 
placing flow symmetry conditions to the lateral side bound- 
aries of the computational domain. 

Results and Discussion 

Test Case A.2 

The A.2 test case was an axisymmetric single-engine 
model with a 0.105m diameter, 0.993m length body mounted 
0.382m from the tunnel wall on a 10 percent thick strut. 
Data were obtained in the Center. The A.2 case consisted 


of a cold jet condition (A.2.1) and a hot jet condition 
(A.2. 2). The jet total temperature for A.2.1 was T tj = 300K 
and 935K for A. 2.2. Both A.2.1 and A.2.2 were tested at 
free stream conditions of M = 0.8, pt ^ = 100000 Pa. Tt x = 
300K with the jet operating at an NPR = 4.8. Experimental 
data obtained by ONERA include streamwise and axial 
flow velocities, turbulent kinetic energy and Reynold s stress 
surveys at several stations downstream of the nozzle exit. 
Only data from A.2.2 is presented in this paper. 

The computational blocking used to simulate the A.2 
model is shown in figure 1. A total of 27 blocks were used 
for the complete configuration. An in-house algebraic grid 
generator was used to generate the single-cell axisymmetric 
grid. The grid included modeling of the free stream in front 
of the model; the model nose, centerhody, and boattail; 
and nozzle internal geometries with high grid density in 
the nozzle trailing edge region. The nozzle trailing edge 
shape was difficult to accurately define due to uncertainty 
in the mating of the internal nozzle and external afterbody 
shells, as well as, some thermal expansion movement due 
to the hot jet. The original surface geometry definition 
used was that lofted by McDonnell-Douglas Aerospace in 
St. Louis, Missouri. A detail of the trailing edge grid is 
shown in figure 2. The modified definition, figure 3, was 
developed later in the study after conversations with other 
program participants. The modified shape appears to be 
more representative of the actual configuration. 

Figure 4 is a comparison of the experimental nozzle boat- 
tail surface static pressure distributions with computed data 
using the two different boattail definitions. Several modifi- 
cations to the boattail definition occurred in the process of 
modifying the trailing edge shape. The transition between 
two sequential conical sections defining the boattail geom- 
etry was sharpened resulting considerably stronger flow ac- 
celeration from about ~1.8 < x/D m < —0.6. The terminal 
boattail angle remained unchanged and as a result the static 
pressure distribution from — 0.6 < x/Dm < 0 were very sim- 
ilar for both trailing edge shapes. The over- stagnation of 
the flow at the nozzle trailing edge is typical of turbulent 
flow calculations using the Jones k Launder k-e model. In 
general, the impact of the trailing edge geometry on pre- 
dicted nozzle boattail static pressure distribution was small, 
especially when compared to the effect of corner sharpness 
between the two conical boattail sections. 

Figures 5 and 6 are flowfield Mach contours in the vicin- 
ity of the nozzle/ boattail for the original and modified trail- 
ing edge geometry configurations, respectively. Figures 7, 8, 
and 9 are radial distributions of several flowfield quantities at 
station (survey plane) x/D m = 0.4840. The vertical dashed 
lines in figures 5 and 6 are the approximate locations of the 
survey plane. There were some small detailed changes in the 
plume which resulted from the t.e. geometry modification. 
The peak Mach number at the plume centerline increased 
from M = 2.58 ( x/D m = 0.42) for the original grid to 
approximately M — 2.84 for the modified grid due to the 
stronger expansion from the inner lip of the modified geom- 
etry. The comparison of computations with axial velocity 
LV-survey data at x/D m = 0.4840, figure 7, shows excel- 
lent overall agreement using either geometry, except near 
the plume centerline. Figure 8 presents a comparison of 
predicted axial normal stresses \u ( \/uoc with experimental 
data. The predicted peak is lower than the data which is 
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Table 1. Test Cases 


Identification 

Free-stream conditions 

Jet flow condition 

Position 

Case 

Horiz. 

tail 

Vert. 

tail 

#of 

engines 

Nozzle 

shape 

Boattail 

separation 

Moo 

Nom. 
?oo, psi 

Nom. 

°R 

NPR 

Nom. 
Tt r °R 

q, deg 

B.1.1 

None 

None 

One 

Axi. 

Yes 

0.901 

4.94 

590 

2.033 

530 

0.010 

B.l. 2 

None 

None 

One 

Axi. 

Yes 

0.901 

4.94 

590 

4.991 

530 

0.009 

B.2.1 

Aft 

Fwd 

One 

Axi. 

Yes 

0.900 

4.94 

590 

2.027 

530 

0.021 

B.2 .2 

Aft 

Fwd 

One 

Axi. 

Yes 

0.900 

4.94 

590 

2.997 

530 

0.021 

B.2.3 

Aft 

Fwd 

One 

Axi. 

Yes 

0.900 

4.94 

590 

5.033 

530 

0.021 

B.3.1 

Mid 

Mid 

Two 

Axi. 

Yes 

0.901 

4.94 

590 

3.407 

530 

-0.004 

B.3.2 

Mid 

Aft 

Two 

Axi. 

Yes 

0.899 

4.94 

590 

3.400 

530 

-0.004 

B.3.3 

Mid 

Fwd 

Two 

Axi. 

Yes 

0.899 

4.94 

590 

3.402 

530 

0.031 


B.1.1 

Cdpn 

Cdpa 

Cdpt 

Cdtt 

Experiment 

-.0015 

0.0291 

0.0276 ± 0.01 

0.0500 ± 0.005 

PAB3D-V10 

-.0209 

0.0260 

0.0051 

0.0389 

PAB3D-Mod.t.e. . . . 

.0006 

0.0281 

0.0287 

0.0654 
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B.1.2 

Cdpn 

Cdpa 

Cdpt 

Cdtt 

Experiment 

-.0096 

0.0296 

0.0191 ± 0.01 

0.0423 ± 0.005 

PAB3D-V10 

-.0330 

0.0253 

-0.0077 

0.0259 


likely caused by the linear k-e turbulence model formulation 
in which there is no separation of the normal stresses. The 
actual axisymmetric jet flow does have some anisotropy in 
the viscous stresses that were not modeled in Version 10 of 
the code. Also, no components of normal stress were pre- 
dicted close to the plume centerline ( z/D m < 0.15). This 
discrepancy is also probably caused by the simplification of 
the turbulence modeling. The thickness of the shear layer 
calcul ated is fairly close to that of the data. The Reynold’s 
stress u'w'/u^, figure 9, is predicted fairly well with only 
minor differences with changes in the t.e. geometry. The 
centerline distributions of the axial velocity, shown in fig- 
ure 10 indicate that the modified t.e. caused a higher gradi- 
ent in axial velocity from 0. < x/D m <0.5. This trend was 
observed previously in the Mach contours of figures 5 and 6. 
The larger internal lip radius of the modified trailing edge 
altered the point of the jet shear flow origin. This can be 
seen by the slight shift in the velocity curves at x/D m = 0. 
This slight streamwise shift in addition to the steeper veloc- 
ity gradients made an accurate comparison of the radial flow 
surveys difficult. 

Langley B Test Cases 

Photographs of the Langley test cases B.2 and B.3 in- 
stalled in the Langley 16 -Foot Transonic Tunnel are shown 
in figures 11 and 12 respectively. Test case B.l is identical 
to test case B.2 with the vertical and horizontal tails re- 
moved. Table 1 identifies the empannage arrangements and 
flow conditions for these three test cases. 

Test Cases B.l and B.2 

Axisymmetric, single-engine afterbody experimental data 
with tails off (B.l) and tails on (B.2) were obtained using 
an air-powered, cold jet simulation system in the 16-Foot. 
Transonic Tunnel at the Langley Research Center (see fig- 
ure 11). All comparisons were made for data at M — 0.9. 


Experimental/ computational comparisons were made for the 
tail-off and staggered-tail configurations NPR = 2 (B.1.1 
and B.2.1) and NPR = 5 (B.l. 2 and B.2.3). The interactive 
program Gridgen was used to generate a vertical half-plane 
symmetric five block grid for the configuration with tails. 
A single-cell wide version of the grid was used for the B.l 
(no tails) configurations. A tabulation of the experimental 
and predicted pressure drag and total drag (pressure + skin 
friction) for B.l is shown in the tables above. These drag 
coefficients have been normalized by the maximum cross- 
sectional area of the body. 

B.1.1 - Figures 13 and 14 show comparisons of theory 
with experimental data for the single jet, axisymmetric ge- 
ometry. In general the surface static pressure distributions 
compared well with the data except for the trailing edge por- 
tion of the nozzle where the code over-predicted the stagna- 
tion at the trailing edge (solid line). This pressure coefficient 
mis-match is due in part to the approximation made in the 
nozzle-boattail geometry which modeled the trailing edge as 
a sharp point rather than the combination blunt, base/open 
cavity of the actual model. This conclusion was reached as 
a result of further calculations performed by a summer re- 
search associate at Langley, Kristina Alexander, after the 
initial study was completed. Several variations on model- 
ing the trailing edge region of the nozzle were attempted. 
The most realistic geometric description involved modeling 
the individual blunt bases of the outer afterbody (0.010 in.) 
and inner nozzle (0.025 in.) with the metric break cavity 
that occurred between the two shells. The total trailing 
edge thickness was approximately 0.101 in. thick. Compu- 
tational results on this more realistic trailing edge geometry 
are shown in figures 13 (dashed line) and 14 (solid sym- 
bols). The lower recovery pressure shown in figure 13 for 
the realistic trailing edge increased the predicted drag on 
the nozzle (see tabulated data above) and improved the com- 
parisons with experimental pressure drag coefficients shown 
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in figure 14. Additionally, the 2-equation linear k-€ turbu- 
lence model formulation of Jones & Launder typically pre- 
dicts static pressure levels in the vicinity of a nozzle trailing 
edge that are too high. Improved turbulence modeling has 
the potential of correcting the trailing edge pressure recovery 
discrepancy further. 

B.2.1 - A comparison between the computed external 
static pressure coefficients and experimental data at several 
angular stations for the B.2.1 (staggered tails, NPR — 
2.027) case is shown in figures 15, 16, and 17. The row 
<t> = 0 is along the top centerline of the model in line with the 
vertical tail. The first C p rise (around x/f = 0.64) is caused 
by pressure stagnation at the leading edge of the vertical tail. 
Further downstream, there is a region where the pressure 
distribution was affected by the flow over the vertical and 
the horizontal tail surfaces as well as the afterbody and 
nozzle surfaces. The computed C p values were in general 
good agreement with the measured data. However, the local 
minimum near x/l = 0.9 was not captured in the numerical 
simulation. The computed pressure recovery downstream 
of x/l = 0.9 followed the measured data closely until near 
the end of the nozzle surface. The computed pressure 
coefficient was approximately 10 percent higher than the 
experimental data at the position of the last pressure orifice. 
This discrepancy is probably caused by the same factors 
as discussed for the the over-prediction noted for the B.1.1 
case. Qualitatively, comparisons at other rows of pressure 
measurement positions were similar or better than those 
shown in figure 10. 

A comparison of the integrated afterbody and nozzle 
pressure drag coefficients with measured values are shown 
in figure 18. The following table provides a tabulation of 
the results. These drag coefficients have been normalized by 
the maximum cross-sectional area of the body. 

\ 


B.2.1 

Cdpn 

Cdpa 

Cdtt 

Experiment . . . 
PAB3D-V10 . . . 

-.0130 

-.0304 

0.0410 

0.0426 

0.1055 ± 0.005 
0.1163 


The computed pressure drag coefficient of the afterbody 
compared well with the measured value. The comparison of 
the pressure drag coefficient on the nozzle external surface 
was less favorable. As one may expect from the over pre- 
diction of the pressure recovery in this region, the pressure 
recovery benefit was also over-predicted. 

Test Case B.3 

The twin-engine afterbody, cold-jet test data were ob- 
tained in the NASA-Langley Research Center 16-Foot Tran- 
sonic Tunnel, (see figure 12). Surface static pressure data 
were measured on the afterbody and tail surfaces. Nozzle 
pressure drag and total aft-end (includes afterbody, noz- 
zle, and tail surfaces) drag are available at M = 0.9 at 
NPR = 3.4. The program Gridgen was used to generate 
the computational meshes. Only the mid-vertical tail (B.3.1) 
and the forward- vertical tail (B.3. 3) configurations (see Ta- 
ble 1) were gridded and approximately 1.5 million grid points 
in 11 blocks were used to describe the model including the 
nose, wing, afterbody, and nozzle internal geometry. 


The computed and measured pressure coefficients along a 
row at the top-center position above the left side engine are 
compared in figure 19. Both the mid-position (B.3.1) and the 
forward position (B.3. 3) vertical tail leading edge positions 
are clearly indicated in the C p distributions (x/l = 0.782 
and 0.727, respectively). In the region between the vertical 
tails, the flow expanded steadily to attain supersonic speeds 
in both cases with C p values as low as —0.56. The flow 
expansion was terminated by a shock ahead of the trailing 
edge of the vertical tail in both cases. An unusual detail at 
x/l = 0.95, in the form of a small dip in C p was captured in 
these flow solutions. This feature is thought to be caused by 
a slope discontinuity between the afterbody and the nozzle. 
Pressure recovery on the nozzle compared well with data. 

The pressure coefficient comparisons along the model top 
surface centerline on the interfairing between the engines 
and nozzles are shown in figure 20. The last data point 
at x/l = 0.95 was near the trailing edge of the interfairing 
between the twin engines. The computed flow expansions 
and shock positions were similar to those shown in figure 19. 
This result indicates that the flow expansion and shock near 
the vertical tail surfaces (and noted for figure 19) carried all 
the way across the entire afterbody upper surface between 
the twin vertical tails. Similar to the B.2.1 case, there are 
several other pressure orifice rows where C p distribution 
comparisons between computation and experimental data 
were made for test case B.3, but are not shown in the current 
paper. The agreements were similar in quality and in some 
cases better than those shown in figures 19 and 20. 

Some detail of the B.3.3 flow solution on the inboard side 
of the nozzle is shown by the computed particle traces in 
figure 21. Although the pressure distribution on the inboard 
side of the nozzle was similar to that over the top of the 
nozzle, a mild flow separation had occurred, according to 
the computations. The computed isobars in the tail region 
of the B.3.3 case are shown in figure 22. The pressure 
gradients at the leading edge of the vertical tail, the flow 
expansion and shock between the vertical tails surfaces, and 
pressure recovery on the nozzle surface are clearly shown in 
this figure. 

Measured and computed nozzle pressure drag and total 
aft-end drag for twin-engine test cases B.3.1 and B.3.3 are 
tabulated in the following table. These drag coefficients have 
been normalized by the wing reference area of the model. 


Conf. 

Cdpn 

Cdtt 

Exp. B.3.1 

-.0020 

0.0084 ± 0.0005 

PAB3D-V10 .... 

-.0020 

0.0106 

Exp.B.3.3 

-.0022 

0.0078 dt 0.0005 

PAB3D-V10 .... 

-.0022 

0.0102 


The measured nozzle pressure drag values were predicted 
well by the computation but the total aft-end drag val- 
ues (includes skin friction predictions) were not well pre- 
dicted. However, it is interesting to note that the predicted 
drag increment between B.3.1 and B.3.3 was well predicted. 
This would indicate that this level of grid density and the 
V.10 simplified k-€ turbulence model is not able to predict 
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absolute drag but may be useful in predicting effect on drag 
of incremental configuration changes. 

Conclusions 

The results of this investigation are summarized as 
follows: 

\ 

1. Jet plume characteristics were matched closely by com- 
putational methods. The effect of trailing edge geometry 
modeling was discernible, but in general did not alter the 
comparisons with data. Limitations in linear turbulence 
models will preclude prediction of any anisotropy in the 
plume turbulence. 

2. External static pressure distributions on single-engine 
and twin-engine configurations were predicted quite well 
matching both level and location of the experimen- 
tal data features. Improvements in the single-engine, 
axisymmetric-nozzle trailing edge pressure recovery was 
obtained by a more accurate grid model of the geometry 
with emphasis in the trailing edge region. Improved tur- 
bulence modeling potentially could further improve the 
predicted pressure recovery. 

3. Compared to experiment, integrated pressure drag pre- 
dictions were typically better than total drag predictions. 
This level of grid resolution and simplified linear turbu- 
lence model may preclude accurate total force predictions. 
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Afterbody detail 




Figure 7.- Comparison of axial velocity with data, 
A.2.2. 



Figure 9.- Comparison of u’w’ Reynold’s stresses 
with data, A.2.2. 



Figure 8.- Comparison of axial normal stresses with 
data, A.2.2. 



Figure 10.- Effect of t.e. geometry on centerline 
axial velocity. 


8 

American Institute of Aeronautics and Astronautics 



Figure 12.- Photograph of B.3 in the Langley 16-Foot Transonic Tunnel. 
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Figure 1 3.-Comparison of experimental and calculated 
pressure coefficient, M=0.90,NPR=2.033. 


Figure 1 5.- Comparison of experimental and calculated 
pressure coefficient, M=0.90,NPR=2.027. 
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Figure 17.- Comparison of experimental and calculated 
pressure coefficient, M=0.90,NPR=2.027. 


Figure 14.-Comparison of experimental and 

calculated pressure drag, B.1 ,M=0.90. 


C p 


Figure 16.- Comparison of experimental and calculated 
pressure coefficient, M=0.90,NPR=2.027. 



NPR 

Figure 18.- Comparison of experimental and 
calculated pressure drag, B.2. 
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DATA, B.3.1, Model CL 



Figure 1 9.- Comparison of experimental and 

calculated pressure coefficients, B.3. 


Figure 20.-Comparison of experimental and 

calculated pressure coefficients, B.3. 


Afterbody detail 
Tails omitted 
M=0.9 
NPR=3.4 



Figure 21 .-Particle trace on nozzle external surface, B.3.3. 
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